Superconductivity from doping boron icosahedra. 
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We propose a new route to achieve the superconducting state in boron-rich solids, the hole doping 
of Bi2 icosahedra. For this purpose we consider a prototype metallic phase of B13C2. We show that 
in this compound the boron icosahedral units are mainly responsible for the large phonon frequencies 
logarithmic average, (65.8) meV, and the moderate electron-phonon coupling A = 0.81. We suggest 
that this high Tc could be a general feature of hole doped boron icosahedral solids. Moreover our 
calculated moderate value of A excludes the formation of bipolarons localized on the icosahedral 
length scale as suggested by previous authors. 
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I. INTRODUCTION 

Low atomic number elements have been intensively in- 
vestigated in the attempt of finding electron phonon me- 
diated superconductors with high critical temperatures 
{Tc). Boron rich solids are eminent examples. The high 
energy phonon frequencies of metallic boron layers in the 
structure of magnesium diboride^ are mainly responsi- 
ble for the 39 K Tc- Intercalation of these boron lay- 
ers with lithium and boron substitution with carbon are 
currently under study^ in the quest of even higher Tc- 
Elemental boron becomes metallic (in a non-icosahedral 
structure'^''*) and superconducting at 160 Gpa, with a Tc 
which increases up to 11.2 K at 250 Gpa^. A 23 K Tc has 
recently been discovered in intermetallic, yttrium palla- 
dium boron carbides'^. 

High temperature superconductivity has been found in 
alkali intercalated fuUerene (Cgo) and theoretically sug- 
gested in other intercalated carbon polyhedrons, €20®, 
C2s^ or Cae-^. Cgo is a band insulator under normal con- 
ditions but the intercalation with alkali atoms does gen- 
erate a metallic state and a consequent superconducting 
state with Tc up to AQK^. Superconductivity is mainly 
sustained by the high frequency intramolecular phonons 
so that most of the physical properties of alkali doped 
fuUerenes can be understood from the solid Ceo electronic 
and phonon structures. In particular it is seen that the 
small radius of the molecule substantially increases the 
electron-phonon coupling with respect to the case of the 
unrolled graphite layer—. The study of solids made of 
light atoms is then interesting since they can satisfy the 
two important requirements of having large phonon fre- 
quencies and fairly high electron-phonon coupling. 

In this work we propose a new route to achieve the 
superconducting state in boron-rich solids, the hole dop- 
ing of B12 icosahedra. For this purpose we consider a 
prototype metallic phase of B13C2 (sec. ^1 for which 




FIG. 1: Atomic structure of B12C3 and of B13C2. The black 
atoms are the polar sites, bonded to neighboring icosahedra. 
The white atoms are the equatorial sites. The gray atoms 
form the chain. In B12C3 the atoms in the chain are CBC 
and the icosahedra are BuC with the carbon placed in a polar 
site. In B13C2 the carbon atom in the polar site is replaced 
with a boron atom. 

we predict a Tc comparable to that of MgB2 . We show 
(sec. IIII|) that in this compound the boron icosahedral 
units are the main responsible for the predicted large Tc. 
Thus this high Tc is a general feature of hole doped boron 
icosahedral solids. 



II. CRISTAL STRUCTURE 

We consider the hole doping of boron carbide B12C3 
(or B4C), a wide gap band- insulator. Its crystal structure 
consists of an arrangement of BuC distorted icosahedra 
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FIG. 2: DFT band structure of B13C2. The energy (eV) is 
referred to the Fermi energy level. The convention for high 
symmetry points in the Brillouin zone is from ref. 33. 



on the site of a rhombohedral lattice and of a linear C- 
B-C atom chaioiiii^iiSJiiiS-iSiii. The periodic unit cell 
contains 15 atoms and is illustrated in fig. ^ Hole doped 
boron icosahedra can be naturally obtained replacing in 
B12C3 a carbon atom by a boron, giving B13C2. Due 
to the carbon substitution there is one hole per unit cell 
and, B13C2 being non magnetic, band theory predicts a 
metallic behaviour. At zero temperature and ambient 
pressure however, B13C2 is a semiconductor—. The in- 
sulating character might be related to (i) the presence of 
structural defectsi^iSS (U) Mott polaronic featuresSi. In 
the first case a metallic state might be achieved simply 
by producing high quality samples. In the second case 
it might be generated by the application of hydrostatic 
pressure in order to increase the hopping between the 
icosahedral units. The pressure necessary to metallise 
B13C2 would be lower than that necessary to metallise 
B12C3 or a—B. Indeed, the Mott-polaron gap in B13C2 
is expected to be much smaller than the several eV gaps 
of the band insulators B12C3 and a—B. 

The two most probable structures of B13C2 are 
BiiC(BBC), i.e. BnC icosahedra linked by BBC 
chains, or Bi2(CBC), i.e. boron icosahedra linked by 
CBC chains. From the theoretical point of view. Den- 
sity Functional Theory (DFT) calculations^^ identified 
the Bi2(CBC) structure as the most stable one with 
a 2.09 eV/cell larger binding energy with respect to 
BiiC(BBC). Experimentallyii the NMR spectrum 
of B13C2 confirms the DFT calculations. Indeed chang- 
ing composition from B12C3 to B13C2 the NMR peak 
associated to the C atom in the icosahedron disappears. 
Therefore, in this paper, we use the Bi2(CBC) structure 
and we refer to it as the B13C2 crystal structure. In other 
works it has been suggested that the substituted boron 
is on the chain^^. We did not consider this proposed 
structure since we do not expect that the conclusions 
presented in this work are substantially affected by the 
location of the additional B atom. 



III. CALCULATIONS AND RESULTS 

In B13C2 a metallic phase could be achieved either by 
obtaining clean samples or by applying a small pressure. 
In DFT, already at room pressure, the ground state is 
metallic. Therefore we use DFT to study the super- 
conducting properties of the hypothetical metallic state 
of B13C2. Electronic structure calculations^i and geo- 
metrical optimization are performed using DFT in the 
local density approximation. We use norm conserving 
pseudopotentials'^^ with a s non local part. The wave- 
functions are expanded in plane waves using a 40 Ry cut- 
off. For the electronic structure we sample the Brillouin 
zone (BZ) using a 4'^ Monkhorst-Pack grid (10 k-points 
in the irreducible BZ wedge) and first order Hermite- 
Gaussian smearing of 0.03 Ry. From geometrical opti- 
mization we obtain the values of a = 9.686 (a.u.) and 
a — 66.05 (deg.) for the cell parameters of the rhombo- 
hedral unit cell, very close to the experimental valuesii 
(a = 9.823 a.u. and a = 65.62 deg. ). 

The DFT band structure of B13C2 is shown in fig. 
121 and is in good agreement with ref. It displays 

a metallic behaviour. The Fermi level is close to the 
top of the valence band and it is crossed by several 
bands. We calculated the electronic density of states 
(see fig. 0)) of B13C2, using a mesh of Nk = 14'^ in- 
equivalent k— points. The mesh is generated by tak- 
ing the mesh centered at the F point and shifting it by 
a random vector. The density of states at the Fermi 
level is A''(0) = 3.6 states/eV/unit cell. We decompose 
the A''(0) in icosahedral and chain density of states, by 
projecting the ab initio B13C2 wavefunction on the ba- 
sis formed by the respective atomic pseudo- wavefunction 
(Lowdin population)^. At the Fermi level the icosahe- 
dral states are responsible for 88% of the total density of 
states {Nico{0) = 3.2 ). In particular the boron atoms in 
the polar sites have the largest contribution to Nico{0), 
namely Npoiar{0) = 2.3 and the contribution due to the 
equatorial sites is smaller, Neq{0) = 0.9. Thus, most of 
the electrons involved in conduction processes resides in 
icosahedral states and only a small part in chain states, 

iNcha^n{0) = 0.4). 

We compute the harmonic phonon frequencies ujq^i, us- 
ing Density Functional Perturbation Theory in the linear 
response24. We use a iV = 4^^ Monkhorst-Pack q— points 
mesh, q being the phonon wavevector. The total phonon 
density of states, F{uj), together with the phonon density 
of states restricted to the icosahedral and chain phonon 
modes are shown in fig. ^ The large number of peaks 
in F(lu) is determined by the large number of phonon 
modes present in the system. 

The icosahedral phonon modes are responsible for most 
of the weight between and 140 meV. The vibrations of 
the atoms in the chain explain the high energy feature 
at 193 meV (both B and C vibrations) and part of the 
feature at 129 meV (C vibrations). The structures in the 
chain restricted F(uj) between 25 and 55 meV are mainly 
due to B vibrations, while the remaining weight between 
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FIG. 3: DFT electronic density of states of B13C2 compared 
to the icosahedral atoms and chain atoms projected density 
of states. The energy (eV) is referred to the Fermi energy 
level. 



60 and 106 mcV is due to C vibrations. 

The electron-phonon interaction for a phonon mode v 
with momentum q can be written as: 



A 



u:^N{Q)Nk 



l5'kn,k+qmP'5(£kr,)(5(£k+qm) (1) 



k,n,m 



where the sum is carried out over the BZ, and ekn are 
the energy bands measured with respect to the Fermi 
level at point k. The matrix element is g'^^ k+qm ~ 
()in\5V / 5uc^v\\i + cim) / yj2uj^ , where Uqy is the ampli- 
tude of the displacement of the phonon v of wavevec- 
tor q, V is the Kohn-Sham potential and iV(0) = 3.6 
states/eV/unit cell is the electronic DOS at the Fermi 
level. The electron-phonon coupling A is then calculated 
as an average over the N q— points mesh and over all the 
modes, A = Eq^Aq,/iV = 0.81. 

The modes responsible for superconductivity can be 
identified from the Eliashberg function o^FIlo) 



qi/ 



(2) 



The Eliashberg function is depicted in fig. |31 Most of the 
contribution to A comes from the region from 60 to 105 
meV, mainly related to the icosahedral phonon modes. 

The critical superconducting temperature is estimated 
from the calculated phonon frequencies and electron- 
phonon coupling using the McMillan formulaSi: 



M 

1.2 



cxp 



1.04(1 + A) 
A-^*(1 + 0.62A) 



(3) 



where /x* is the screened Coulomb pseudopotential which 
takes into account the Coulomb repulsion between the 
electrons dressed by retardation effects due to the 
phonons and {lo) = 65.8 meV is the phonon frequen- 
cies logarithmic average. The calculated values of Tc for 



FIG. 4: Phonon density of states, F{u), of B13C2. The total 
density of states (solid line) compared to the phonon density 
of states of the icosahedral modes (dashed line) and of the 
chain modes (dotted line). Most of the weight between and 
140 meV is due to icosahedral modes. 
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FIG. 5: Eliashberg function (solid line) and average electron- 
phonon coupling (dashed line) of Bi3C2.The icosahedral 
phonon modes ranging between 60 and 105 meV give the 
largest contribution to A. 



B13C2 as a function of /x* are illustrated in table The 
critical temperature for metallic B13C2 is comparable to 
the one obtained for MgB2 and, using the McMillan for- 
mula, ranges between 15.8 k and 36.7 K. 

Our calculation of the electron-phonon coupling gives 
also new insight on the possible occurrence of a bipo- 
laronic insulating state in B13C2, as suggested by other 
authors?". In icosahedral boron compounds the bipo- 
laron could be localized on two length scales: the icosa- 
hedron length scale or a smaller one, of the order of the 
bond length. If the bipolaron were localized on the icosa- 
hedron length scale (as most of the bipolaronic litera- 
ture seems to suggest^^) the bipolaronic distortion could 
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Material 


A'* 


(tj) (meV) 


A 


Tc (K) 


B13C2 


0.1 


65.8 


0.81 


36.7 




0.14 


65.8 


0.81 


27.6 




0.2 


65.8 


0.81 


15.8 


MgBa 


0.14 


62.0 


0.87 


30.7 



TABLE I: Critical temperatures of B13C2 and MgB2. Pre- 
dicted critical temperatures as a function of the screened 
Coulomb pseudopotential (/i*). The critical temperature is 
estimated using the McMillan formula (eq. EJ.The results are 
compared with MgB2 (from ref.'^^). 

be seen as a small perturbation to the metallic state. 
In this case the bipolaron could be described within a 
model considering a linear electron-phonon coupling per- 
turbation to the harmonic Hamiltonian in the metallic 
phase. Such a model is the one used in the present paper 
to study the superconducting properties. Our calculated 
value of A is too small2^ to justify a bipolaronic insulating 
state for B13C2. On the contrary our calculation cannot 
exclude the occurrence of a bipolaron localized on the 
bond-length scale. In this scenario the bipolaron would 
involve a substantial deformation of our crystal structure 
(e.g. the breaking or formation of chemical bonds) . Such 
a strong deformation could not be seen as a weak pertur- 
bation to the harmonic Hamiltonian in the metallic state 
and could not be described by the approach used in the 
present paper. 



B13C2, which is formed from B12 icosahedral units with 
one hole per icosahedra. Using ab-initio calculations we 
have determined its normal state properties, finding a 
moderate electron-phonon coupling and a large phonon 
logarithmic average to the phonon frequencies. We have 
demonstrated that both properties are connected to the 
B12 building blocks. Indeed the local density of state at 
the Fermi level and the phonon modes strongly coupled 
with electrons are localized on the icosahedra. As a con- 
sequence our findings are not restricted to B13C2 com- 
pounds but can be applied to other metallic compounds 
composed by B rich icosahedra. For example another 
possibility to achieve a metallic state is the substitution 
of P with Si in B12P2 or of As with Si in B12AS2. B12P2 
and B12AS2 are band insulators composed by B12 icosa- 
hedra and 2-atom P2 or As2 chains^'^. The substitution 
of pentavalent atoms like As and P with a tetravalent Si 
introduces a hole in the system. Bi2P2-2:Si2; wafer resis- 
tivity measurements show a low electrical conductivitj^. 
If such conduction were related to bipolaron formation 
then a metallic state could be achieved applying pressure. 
An advantage of the 2-atom chain systems (Bi2P2-a;Sia; 
or Bi2As2-a;Sia:) respect to the 3-atom chain systems 
(B13C2) is the lack in the formers of an internal soft 
degree of freedom, the chain bending. Probably this fea- 
ture makes the 2-atom chain structures more stable under 
pressure. 



V. ACKNOWLEDGEMENTS 



IV. CONCLUSIONS 

In this work we have studied the possible occurrence 
of superconductivity from hole doping boron icosahedra. 
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